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In the Euclidean plane: given three non-colinear points (vertices of the "triangle of 
reference") any point is uniquely determined by its distances to those vertices. These are 
called the "distance coordinates" of the point. The main result of our first paper was to 
determine which vectors of three non-negative numbers could be distance coordinates for 
the given reference triangle. In this paper we put that result, and others, into matrix form. 
This leads to generalizations, and to the effect of a change of the reference triangle on the 
distance coordinates and the formulas in which they are involved. 
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In [1] * we defined distance coordinates for a point in the Euclidean plane with respect to a 
triangle of reference. We exhibited formulas relating these coordinates to the area coordinates 
of the point. 

In this paper we put those formulas into matrix form, note how they are effected by a 
change in the triangle of reference, and generalize them considerably. The results are listed in 
section III, labeled (R.l) through (R.30). 

Sections I and II state our definitions and notations, and contain equations (1) through 
(27). Sections IV and V contain our proofs and equations (28) through (63). 

I. Matrix Definitions 

Unless stated otherwise capital Roman letters will denote matrices of order three. For 
example I is the identity matrix, and J is the matrix with 1 in each position. Small Greek 
letters will denote one by three vectors. For example 

e = (l 1 1). 

Given a matrix X=(x i; ) its transpose will be denoted by 'X T =(x ji ). Its determinant is 
det X. Its adjoint is 



(1) X A = 



^^22^33 £23^32 ^13^32 ^33^12 ^12^23 ^22^13^ 



#23^31 #33^21 ^11^33 # 13^31 #21#13 #11^23 

k^32#21 #22^31 ^31^12 ^11^32 X\\X 22 X\ 2 X 2 iJ 



Of course 

(2) XX A =(detX)I=X A X. 
The sum of the elements in X A is 

(3) h(X)=x n X22+XnXz3+x 2 2X33— x u x 2 i — x lz x n — x 2Z x Z2 — ^11(^23+^32)— ^22(^13+^31)— ^33(^12+^21) 

+ ^12^23 + ^13^32 + ^21^13 + ^23 + ^31^31^12 + ^32^21 

Note that (X A ) T =(X T ) A and that h(X T )=h(X). 
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Given a vector ti=(m,i m 2 m 3 ) we define the matrix 

f m 3 — m 2 ' N 



s,= 



-m 3 mi 
^7712 —mi „ 



It has the properties 



det S M =0, mS m =0 7 S m t =-S„ 
S m a =mV 



and 

Also 

(4) ^S M = — juS„ 

for any vector v. The special notation S is used for S e 

r 1 -r* 



-1 1 
1 -1 Oj 



We have 

(5) 



^A /£, 



=J. 



Throughout this paper we shall assume the following results which we shall prove in the 
last section. 

Given a matrix X and vectors /x, v we have 

(6) det (X+ M T .)=det X+yXV 
and 

(7) (X+ M T ,) A =X A +S,X T S M T . 
Given matrices X, Y, Z, with A(X)?^0, such that 



(8) 

then 

(9) 

If h(Y) ^0 then 

(10) 

Finally for any matrix X : 

en) 

or, equivalently, 
(12) 



/t(X)Z=JX A +YSX T S T , 
det Z=A(Y)/A(X). 

A(Y)Z- 1 =JY A +XSY T S T . 

XSX T S T =A(X)I-JX A , 

SX T S T X=A(X)I-X A J. 
62 



II. Geometric Definitions 



Let t be a triangle in the Euclidean plane with area A,>0 and vertices PuP 2 ,P3- With re- 
pect to t (the triangle of reference) a point p has unique area and distance coordinates: 

P=(bi,b 2 ,b 3 )=[d u d2,<k] 

respectively. The b t are normalized (i.e. &i+& 2 +&3=1) barycentric coordinates, and 

d t =\p— Pi\ i=l, 2, 3 

is the distance from p to the vertex p t . Define the vectors 

Pp/t= (616263) 
and 

dpn=(di 2 d2 2 <k 2 ). 

For example for the vertices of t we have 

fttf=(l 0) 

(13) fr*=(0 1 0) 

ftw*=(0 1) 

and 

8 Pl /t=(0 a 3 2 a 2 2 ) 

(14) 8 P2/ t=(a 3 2 a, 2 ) 

d P3 /t=(fl 2 2 a, 2 0) 

where a { is the length of the side of t opposite the vertex p u i=l } 2, 3. 

Let t\ be a triangle with area A^ and vertices q u q 2 , 2a- Define the matrices 



"Btdt= 



An equivalent definition for D tl / t is 

(15) Vh!t=(\qi-Pi\ 2 ), 
so that 

(16) ~D m =T>? llt . 
For example, equations (13) tell us that 

and equations (14) tell us that 

r a 3 2 a 2 2 ^ 

(17) D„,= a 3 2 a, 2 

^a 2 2 a? J 



l%/< "I 




V 


£W, 


D^= 


<W 


<Pq2/t > 




J>Qzftj 
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We note that 

(18) 
that 

(19) 

that 
(20) 
where 

and that 

(21) 



det D t/t =2a 1 2 a 2 W, 



D A — 



r —a 1 i a 2 a 2 2 a^a^ 
a^a 2 2 — a 2 4 a 2 2 a s 2 



^&l #3 d 2 &3 &3 ^ 

eD? /t = (<ii 2 Ci a 2 2 c 2 dz%) 



c i =a 1 2 +a 2 2 +a 3 2 —2a-' 



i=l,2,3 



/.(D^eD^e^lGA, 2 . 



It is also important to note that 



(22) 



SD^S T = [c z 



-2d! 2 



c 3 



c 2 

Ci 

-2a, 2 



For any point p define g plt to be the negative of the power of p with respect to the circum- 
circle of t. That is, g p/t is the square of the circumradius of t minus the square of the distance 
from p to the circumcenter of t. Define the vector 



(23) 



yti/t=(gqiitgq2/tgqz/t)- 



Note that g p /t=0 if and only if p is on the circumcircle of t. Since each vertex of t is on its 
circumcircle, it follows that y tjt =Q. 

To the results listed in the next section should be added the following which are well known 
(see [2], pages 218-219). 

(24) |detB^|=A^. 
Thus 'B h /t is non-singular if and only if A^>0. In that case 

(25) B^ = B^, 
and 

(26) Bp/*, = Bp/t Btfti 

for any point p. Moreover, if t 2 is any triangle with A; 2 >0 then 

(27) B tl /t Bf/fc=Bfc# 2 . 

Finally, if /3 is any real one by three vector such that /5e T =l then there is a point p such 
that p=p p/t . 
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III. List of Results 

If 5 is any real, non-negative one by three vector, then there is a point p such that 8 = 8 p/t if 
and only if 

(R.l) det(D t/t -e T 8-8 T e) = 0. 

If A^=0 then ~D tl /t is singular (as is B^). But 

(R.2) if A^>0 then D tl /t is singular if and only if the circumcircles of t and t x are orthogonal; 

i.e. the circles intersect and the respective radii to a point of intersection are perpendicular. 
For any two points p, q we have 

(R.3) \v-q\ 2 ={*vi-hit)0lit> 

Thus 

(R.4) («w«-*f/i) (P P/t +P Q/t ) T =o. 

Also 

(R-5) P P /t£>t l /tP 1 Q/t =\p— ql'+gpih+og/t A/,>o. 

In some cases we may desire a formula to be independent of the dimensions of t, or inde- 
pendent of the distance or area coordinates. The following three equations show that this can 
be done for g p/t : 

(R.6) g v it=h/S T pit 



(R.7) 9v/t=oPp/t'Dt/tfi 



2 



J p/t 



(R.8) gpft^ilfo* 2 )- 1 (det D t/ -5 P/t D$ /t € T ). 

We can replace 8 P/t in a formula by using 

(R-9) Spit=$p/tE>tit—(Jpitt- 

And we can replace $ vll in a formula by using 

(R.io) fi p/t =(mA t 2 y^i>tf t +8 p/t SB t/t s T ). 

If we wish to change our triangle of reference from t to t u it turns out that we require 
A(D^)t^C). The following formulas show that this is equivalent to A^>0: 

(R.ll) A(D^) = 16A, 2 detB^ 

(R.12) \h(D tl/t )\ = lQA t A h . 

The first three "change of basis' ' formulas are 

(R.13) 8p/ti=Pp/tDuit—g p /te 

(R.14) h(D tllt )0 p/tl =eD£ /t +8 p/t S D? l/t S T A tl >0 

(R.15) h(D tl/ t)g m =det D^-^,D^e T A fl >0. 
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Note that the last two equations depend only on distance coordinates. (For reliance on area co- 
ordinates alone equations (25) and (26) tell us that P P /t l = Ppit B^ t if A^>0.) 

These three equations can be rewritten 

(R.16) «pfc=fttf (D h it-e T d p/t ) T 

(R.17) h(D tl it) p P / h =e(D h/t -e T 8 p j t ) A A*>0 

(R.18) h(D m ) g p/h =det (D tl/t -e T 5 p/ t) A*>0. 

We note that h(D tl it—e T 5 P it)=h(D tl it). Also 

(R. 19) det (D tl it-e T d p/ t-d p / h e)=0. 

Other " change of basis' ' formulas are 

(R. 20) 8 p /t-8 P it=Pp/t (D tllt -D tlt ) T 

(R. 21) T> tllt -D tllt =B tllt (D tdt -T> m ) T 

(R.22) D^=B^D^- 7 ^€ 

and 

(R. 23) h(D tl/t ) B^=J D&+D# S D^ S T . 

Recall that D^=D^ (equation (16)), and that B t / tl =B^ft if A^>0 (equation (25)). 

The similarities between many of these equations is due to the fact that they are special 
cases of the following formulas. 

Let tu t 2 , tz, ti be any triangles. Then 

(R.24) B m (D, 2/t -Dt 2/ t i ) = 1 D h/ t-T> tl/ti A 6 >0 

(R.25) h(D m )B tl/t2 =J Df 2jh +D m S D^ 3 S T AtA h >0 

(R.26) h(D m ) (B tl/t -B t2lh ) = (D tl/t -B t2/ti )S D^S T AtAu>0. 

Also 

(R.27) B U /t 2 ~Dt2!t i =~D tl jt i J ry? l/t2 e A* 2 >0 

( R - 2 8) IW /(3 =7^-7^ A tAl >0 

(R-29) B './< (l =-%T /(2 ^A>o 

and 

(R.30) B^ 3 BJ ( =D^+7^e+6^ A^ 3 >0. 

For example (R.5) is a special case of (R.30) as follows: 

Reduce triangle ^ to a point p and triangle t A to a point q. Then 

B^ 2 ==e T 0p/; 2 , B^ 3 =€ T ft/^3, Dfc/* 4 =|#— g| 2 J, y tlit2 =g p/t2 e and yt 4 /h=gpjh€. All terms reduce to mul- 
tiples of J=e T e whose equality is (R.5) with t 2 -^ti and 2 3 — >£. 
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IV. Proofs 

We will refer to several equations in [1] and will use the notation "(12')" for "(12) in [1]", 
or "equation (12) in [1]." 

Equation (12') can be written 

3 3 

32A, 2 /"(zi, x 2 , x 3 )=2 XI <*>k 2 c k x k — 2a 1 2 a 2 2 a 3 2 — 2 X) d 2 kX k 2 +2c l x 2 :r z +2c 2 x 1 x 3 +2c 3 x 1 x 2 ' 

k=l k=l 

Let 5=(xiX 2 x 3 ) and look at equations (20), (18), and (22). Then we have 

(28) 32A,70z 1 , x 2 , x 3 )=2eDf /t 8 T -det D t/t +8SD t/t S T 8 T 

We shall prove that the r.h.s. is— det (D t/t — e T 8— 5 T e). Since the Theorem in [1] states that 
(assuming <5>0) there is a point p such that 8 = 8 p/t if and only if /(#i, x 2 , x 3 ) = 0, we will have 
proved (R.l). 

From equation (6) we have 

(29) det (D tJ -e T 8-d T e)=det (D t/t -e T 8)-e(D t/t -e T 8) A 8 T 

= detD t/ -8Vt /t e T -e(D t/ -e T 8) A e T . 



From (7) we have 



e(D t /t—e 8) =eDt/t— eS«D^S 

= eD^+5SD, /f S T 



because D t /t=T>tft and eSa= — <5S (see (4)). Substituting into (29) yields 

det (D t/ t-e T 8-8 T e)= det D^-5D^€ T -eD,> T -5SD^S T 5 T . 

Since 6D^€ T =(eD^6 T ) T , the two scalars are equal. Comparison with the r.h.s. of (28) proves 
our result and thus (R. 1). We shall prove (R. 2) at the end of this section. 

Given points/?, q let I3 p/t = (b l b 2 b 3 ) and fi g/t = (b [ b' 2 b' 3 ). 
From (28') we have 



(30) 




2\p-q\*=T;c k (h-K 


From (33') we have 






(31) 




3 

2fif,/i=Sc*(6*— 6J). 

fc=i 


Thus 






(32) 


2|? 


3 

>-q\ 2 +2gv/t+2g g /t=m c k (b k - 



fe=l 

Using c k =a ] 2 +a 2 2 +a 3 2 — 2a k 2 , fc=l, 2, 3, this becomes 

(33) \p—(l\ 2 +g P / t +g Q /t=ai 2 (b 2 bs+hK)+a 2 2 (bX+hb[) + 
The r.h.s is just P P /iDt/t0q/t so that 

(34) \p-q.\ 2 =P P /tDt/tfi/t—gp/—g q /t> 
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We proceed to show that l3 p/t D t/t p T Qft =g p/t +8 p / t l3 T Q/t . We have /3 p/t D t/t = (b 2 a 3 2 +b 3 a 2 2 ). 

From (43') this becomes ff p /^Dt/t=W+g P f t ) or 

(35) Pp/tD t /t=&p/t+gv/tt- 

Since effi /t = l we have 

(36) 0p/tD t /t(tf/t = &p/t&q/t-\~gp/t- 

Substituting this into (34) we get 

(37) \p-i\ 2 =8p/tfi/t-g Q /t. 

From (37 ; ) we get 



T 



(38) Qpit — hitfipit- 

Then (37) becomes \p— q\ 2 = 8 p /t Pq/t—8q/t Pqft, and (R. 3) is proved. Interchanging p and q yields 
\p—q\ 2 =(d q jt—d p /t) Pp/t. Equating the two r.h.s. proves (R. 4). And (R. 5) was proved by the 
remarks at the end of the last section (once (R. 30) is proved). 

We have already proved (R. 6) with (38). Setting q=p in (36) and using (38) proves (R. 7). 
Multiply (35) on the right by D^ e T to get det T) t /t = d P it Dfjt e T + h(T>tit) g P /t. Using (21) com- 
pletes the proof of (R. 8). 

We have already proved (R. 9) with (35). Equations (16') can be written 

(39) 16A, 2 b x = a x 2 d - 2a, 2 d, 2 + c 3 d 2 2 + c 2 d 3 2 etc. 
Equations (20) and (22) trans 1 ate this into 

(40) 16A r 2 ft,,, - e D&+ h p! t SDtit S T 
which is (R. 10). 

Setting p=q h q 2 , q 3 in (40) yields 

(41) 16A, 2 B ult = JD&+Dm SD# S t . 

Since 16A 2 =h (Dt/t) and D^ = D^ we can apply the results following equation (8) withX=D^, 
Y=D^andZ=B^. Equation (9) proves (R. 11) and equation (10) proves (R. 23) since B tl /t= 
B'tl/t if A^>0. Result (R. 12) follows from (R. 11) and (24). 

Interchange p and q in equation (37) and rewrite it as \p— q\ 2 — Pp/t&H/t—Qp/t- Letting q=q\, 
q 2 , q 3 yields (R. 13). Multiply (R. 23) on the left by ft,,,. Since ft,/3*/*i=ft>/*i, ft/ t J=e, p p/t D t/t = 
8 p/t +g p/t €, and eS=0 we get (R. 14). Result (R. 15) follows from (R. 13) and (R. 14) with the for- 
mula Qpitx—Ppitfipity In the calculation use S T e T =0 and SD^S T D^==A(D^)I— D^J (from (12).) 

Result (R. 16) is (R. 13) rewritten to conform with the two following results. They in turn 
follow from (R. 14) and (R. 15) by using the same kind of calculations as those surrounding 
equation (29). Multiply (R. 17) on the right by e T , note that p p / tl e T =l and that the r.h.s. is 
h(D h / t — e T 8 p / t ) to show that the latter equals h(Dti/t)- To prove (R. 19) write (R. 13) as 

Pp/t(D tl /t— e T 8 p/t — 8p/t^) T =0. 

Since p p/t cannot be the zero vector (because ft,/,€ T =l) it follows that the matrix T>t l /t—e r 8 p / t 
—8 p /t 1 € cannot be non-singular. 

Result (R. 20) is (R. 13) minus (R. 9). Then p=q h q 2 , q z , in (R. 20) proves (R. 21). 
Similarly (R. 22) comes from (R. 13) by setting p=q h q 2 , <Z3, and recalling the definition (23). 
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Let ti=t z in (R. 13) and then setp=q u q 2) g 8 , to get 

with t = t 2 (thus requiring A^ 2 >0) and D^ 2 =D^ 3 we get (R. 27). In that equation set t z =t i 
and subtract from the original to get (R. 24). Let ti=t 2 in (R. 14) and then p=q u q 2 , ft- This 
yields (R. 25) with t = t z . 



Result (R. 26) comes from the following observation: Let X=D^ 2 in equation (11) to get 

JD^=A(D W(3 )I-D (2/i3 SD ( T /(3 S T . 
Substituting this into (R. 25) yields 

(42) A(D^ 3 )(B,/, 2 -I) = (D, l/ , 3 -D^ 3 )SD^ 3 S T . 
Subtract this from the same expression with ti^t^ to get 

h(Dt 2 /h) (Bu/t 2 — Rh/ti) = (Dtt/ti—'Dtiiti) SD^ 3 S T . 
The cyclic permutation (12 3 4) transforms this into (R. 26). 
Use the cyclic permutation (1 2 3) on (R. 27) to get 

B^Dij/p D* 2 /£i + yt 2 /t 3 e A* 3 ^> . 
Assume A <2 >0 and multiply on the left by B^ 2 . Using (27) we get 

IWfcD* 3 /h = B u /t 2 T>t 2 /t l + B^/fcYh/f 3 e . 

Now &t i /t£>tz/ti='Dti/ti+yu/he and B U /t 2 Dt 2 /t l ='Du/u+yu/t 2 e (using (R. 27) with the proper permuta- 
tions). Thus 

lu/h e = 7u/t 2 € + ttujt 2 yt 2 /h € • 

With t 4 =t l the first column of this matrix equation is (R. 28). In (R. 28) set t z =U and recall 
that y*,/fc=0 to get (R. 29). 

Assume A <3 >0 and multiply (R. 27) on the right by B^ a to get 

(43) B tl /t 2 D t2 /t 3 ~Bt i it 3 =T)t 1 /tiBujti-\-ytiit 2 e A* 2 Aj 3 >0 

Using (R. 27) we have B u/ti Dh/t l = 1 Du/h+yu/he. Transposed this is D^ 3 B^ 3 =D^+e T Y^3- Substitu- 
tion into (43) yields (R. 30). 

We finish this section by proving (R.2) in the course of which we will use capital Roman 
letters which do not denote matrices. Specifically R, and O t will denote the circumradius and 
the circumcenter respectively of the triangle t. The distance from O t to each vertex of t is R, 
so that 

(44) d 0t/ t=R t 2 e. 
Letting p — O t in (R.15) yields 

det D tl/t =h(D tl/t ) (<7o,/<+Rr 2 ) ^>0. 
By definition, g p/h =U tl 2 — \p— O tl \ 2 . Thus 

(45) det D tllt =h(B tl/t ) (R^+R^-IO -O^l 2 ) A„>0. 
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Since ^R t 2 -\-'Rt 1 2 =\O t —Ot l \ 2 if and only if the respective circumcircles are orthogonal, we have 
proved (R.2). 

V. Matrix Proofs 

Let ~X.= (x tJ ) be a matrix of order three, let v=(nin 2 n 3 ), and let 

a = v~K A =(s 1 S 2 S 3 ). 

Since A(Q)=eQ At e T , e T e = J, and S e = S, we see that letting co = /x = e in equations (60), 
(61), and (63) yield equations (8), (9), and (10) respectively (with Y=P, X=Q). Similarly, 
letting fx = v = e in equations (51) and (52) yield equations (11) and (12; resp ctively. Finally, 
equation (55) is equation (6), and equation (54) is equation (7). 

By inspection of (1) we have 

(46) s 1 =X2i(n 3 x 3 2—n2X^)+X22(n 1 x^—n 3 x 31 )+X23(n 2 x zl —n 1 x d 2) 

= —x 31 (n 3 X22—n2X2z)—x 3 2(n l X2 3 —n 3 X2i)~x d3 (n2X2i—n 1 x 2 2) 

(47) s 2 =x 3l (n 3 x 12 — n 2 x 13 ) +x 32 (n 1 x 13 — n 3 x n ) J rx 33 {n 2 x n ~n l x l2 ) 

= —Xi X (n 3 x 3 2—n2X 3Z )— ^12(^1^33 — ^3^31) — x l3 {n 2 x 3l — ^1^32) 

(48) s d> =x n {n 3 x 22 —n 2 x 23 ) +^12(^1^23— ^3^21) +^13(^2^21— ^1^22) 

=—X2i{n 3 x l 2—n2X lz )—x 22 (:n l x n —n z Xn)—X2 Z {n2Xii—n l x l 2) 
Let 71, 72, 73 denote the rows of XS, T , where 

S„ T 
We have 

7l = (^3^12 — ^2^13 ^1^13 — ^3^11 ^2^11 — 711^12) 

72= (n 3 x 22 —n 2 x 23 n 1 x 23 —n 3 x 2 i ^2^21—^1^22) 

73 = (^3^32 — n 2 X 33 ^1^33 — ^3^31 ^2^31 — ^1^32> 

Let pi, p 2 , p 3 denote the rows of X: 

Pi = \ X il X i2%i3) 7v = l, 2, 3. 

Clearly 

Pi7i T =0 i=l, 2, 3. 

From equations (46), (47), (48) we see that 

Sl z =P273 T = = — P372 T 

S2 = Ps7i T = z — Pi73 T 






—n 3 


71 2 


n 3 





— Tli 


n 2 


Tli 


o 



Thus X(XS, T ) T =( Pi7 /) equals 

That is 

(49) 



S3 = PlT2 T =— P271 T - 





s 2 - 


*3 " 


Si 


-s 2 \ 

o'J 


=S,. 


XS.X 


70 


=s 

V 


XA. 



.'A 



Let /i be a one by three vector. Since 

(50) S,S/=( ( tm t )I-m 1 V 
(a general result), we have 

(51) XS,XV=(^V)I-M T vX A 
Taking transposes and substituting X T for X yields 

(52) S,XVX-(vXV)I-X\ t m. 

Let 

W=X A +S,X T S/ 
Then 

(X+ M ^)W=(det X) I+XS,X T S7+ M T ,X A . 

since XX A =(det X)I and i/S„=0. Using (51) yields 

(53) (X+ M T ,/)W=(detX+vXV)I. 

Since this is an identity in the variables X, n, v (or rather, their elements) we conclude that W 
is a scalar times (X+m t ^) a - In the latter the elements of order two in the x ti form X A , as is 
true of W, so that the scalar is 1. Thus W=(X4- M Y) A : 

(54) (X+m t .') a =X a +S,X t S >1 t , 
and from (53): 

(55) det (X+fv) =det X+„XV. 

Let P, Q be matrices of order three, and let u be a one by three vector. In equations (54) 
and (55) we will set X=PS„Q T S„ T and i>=o>Q A . Since S„ A =//V we have 

(PS„Q T S/) A = (S/) A (Q T ) A S„ A P A 

(56) =m t m(Q a ) t ^ t ^P a 

= (coQV)m T coP a . 
(The scalars m(Q a ) t " t and wQV are equal). From (49) we have S„q/=QS„Q t . Thus 

(57) (PS„Q T S/+ M T U Q A ) A = (coQ V) M T a.P A + (QS„Q t ) (S„QS„ T P T )S/. 
From (51) we have Q T S,QS„ T =(a>QV)I-a J T ( u(Q A ) T . Since S„co T =0, (57) becomes 

(58) (PS tt Q T S/+ M T a,Q A ) A = («QV) (m t «P a +QS w P t S/) . 
In using (54) we recall that det S„=0. Thus 

(59) det (PS w Q T S/+ M T a>Q A ) = ( w Q A )((oQV T )M T c»PV 

= (coQV) 2 coPV. 
Suppose a)Q A ju T is not 0, and define the matrix Z by 

(60) (coQV T )Z=PS„Q T S/+ M T coQ A . 

From (59) we have (coQV) 3 det Z= («Q V) 2 («>P V) 
or 

(61) detZ =((oPV)/(«QV). 
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The l.h.s. of (58) is (coQV) 2 Z A so that 

(62) (coQV)Z A =QSJ> T S, T + M T "P A . 
Assume coPVVO. Then Z A =(det Z)Z" X so that 

(63) (coPV)Z- 1 =QS.P T S7+ M T coP A . 
Comparing with (60) we note that the transformation P «-* Q takes Z into Z _1 . 
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